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Neutron star mergers versus core-collapse supernovae as 
dominant r-process sites in the early Galaxy 
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Abstract. The astrophysical nature of r-process sites is a long standing mystery and many probable sources have 
been suggested in the past, among them lower-mass core-collapse supernovae (in the range 8 — 10 Mq), higher- 
mass core-collapse supernovae (with masses > 20 Mq) and neutron star mergers. In this work, we present a 
detailed inhomogeneous chemical evolution study that considers for the first time neutron star mergers as major 
r-process sources, and compare this scenario to the ones in which core-collapse supernovae act as dominant r- 
process sites. We conclude that, due to the lack of reliable iron and r-process yields as function of progenitor 
mass, it is not possible at present to distinguish between the lower-mass and higher-mass supernovae scenarios 
within the framework of inhomogeneous chemical evolution. However, neutron-star mergers seem to be ruled out 
as the dominant r-process source, since their low rates of occurrence would lead to r-process enrichment that 
is not consistent with observations at very low metallicities. Additionally, the considerable injection of r-process 
material by a single neutron-star merger leads to a scatter in [r-process/Fe] ratios at later times which is much 
too large compared to observations. 
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1. Introduction 

The principal production mechanisms of heavy elements 
beyond the iron-peak are known since the classical paper 
by Burbidge, Burbidge, Fowler & Hoyle (B 2 FH, |P3o7|) : 
A major fraction of these elements is formed by slow 
or rapid capture of neutrons on seed nuclei. These nu- 
cleosynthesis channels are denoted as the s-process, if 
the capture of neutrons occurs on timescales slower than 
the /3-decay timescales of the newly formed nuclei, or as 
the r-process, if neutron capture operates on timescales 
much shorter than the relevant /3-decay timescales. Yet, 
although the physical requirements for the occurrence of 
r-process nucleosynthesis are well understood (see e.g. 
Pfeiffer et al. 2001}, the astrophysical nature of the dom- 
inant r-process site is still unknown. 

Ultra metal-poor halo stars were formed when the 
s-process had little time to make any significant con- 
tribution to the interstellar medium (ISM). The abun- 
dance pattern of neutron capture elements heavier than 
Ba [Z > 56) in these stars matches the scaled solar system 



r-process abundances remarkably well. This suggests that, 
contrary to s-process nucleosynthesis, the synthesis of r- 
process elements started early in galactic evolution and 
that the r-process is robust for elements heavier than Ba, 
i.e. originates from a single astrophysical site or at least 
occurs under well defined physical conditions (e.g. Sneden 
et al. I2000al Westin et al. |2"UUU1 Cowan et al. but 
see also Hill et al. EE 



2J- On the other hand, lighter neu- 
tron capture elements (30 < Z < 56) in ultra metal-poor 
halo stars show significant deviations from the scaled solar 
system r-process abundance curve, indicating the possible 
existence of a second r-process source. The diversity of 
r-process sources was first proposed based on meteoritic 
data, which seem to require two r-process sources oper- 
ating on different timescales, namely ~ 10 7 years for the 
production of Z > 56 nuclei and ~ 10 8 years for Z < 56 
nuclei (Wasserburg et al. .' 1 906|j . In principle, these two 
r-process sites could either be two types (frequent and 
less common) of core-collapse supernova (SN II) events 
(Wasserburg & Pian ^M)! QianEDOOl EOOII) , or a mixture 
of SNe and neutron star mergers (Rosswog et al. | 
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2000}. On the other hand, Cameron (|200ip showed that 
it is possible to reproduce the abundance distribution of 
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both light and heavy neutron capture elements in different 
regions of a single source, namely the accretion-disk and 
jet forming near the proto-neutron star in SN II events. 
However, the evidence that the r-process beyond Ba (at 
least in the range 56 < Z < 72) is robust is very con- 
vincing and we will use this property in the following to 
gain some insight into the enrichment of the Galaxy with 
r-process elements beyond Ba. 

A number of possible astrophysical sites responsible 
for the robust r-process were put forth in the past, most 
of them linked to the violent death of massive stars in su- 
pernova events. Wheeler et al. ( 1998) suggested that, dur- 
ing the prompt explosion of a massive star in the range 
8 — 10 Mq, physical conditions in the innermost mass lay- 
ers are sufficiently extreme for the r-process to work and 
that the amount of r-process matter ejected may be consis- 
tent with observed Galactic r-process abundances. This is 
in agreement with detailed r-process calculations based on 
an artificially induced prompt explosion of an 11 Mq pro- 
genitor (Sumiyoshi et al.|2001). However, there are major 
objections to the prompt explosion mechanism from de- 
tailed studies by Brucnn (1989a, 1989b). Optimistically 
speaking, if prompt explosions of massive stars may oc- 
cur in reality, lower-mass progenitors are the only ones 
that show any hope to do so (e.g. Wheeler et al. [1998 
Sumiyoshi et al. 120011 but see also Liebendorfer et al. 

mnj. 

On the other hand, neutrino-driven winds from nascent 
neutron-stars were also proposed as a promising site for 
r-process nucleosynthesis (e.g. Woosley & Hoffman 119921 
Takahashi et al. IT5M1 Woosley et al. IT3911 Qian & 
Woosley HH2D Thom pson et al. lMiTl Wanajo et al. lMiTl 
Terasawa et al. 2002). Following the delayed explosion of 
a very massive star (> 20 Mq), neutrinos diffusing out of 
the contracting proto-neutron star heat and ablate mate- 
rial from its surface and a neutrino-driven wind develops. 
Such winds may exhibit sufficiently high entropies and/or 
sufficiently short dynamical timescales for the r-process 
to occur. However, r-process yields consistent with ob- 
served r-process abundances in stars may be obtained in 
this scenario only for extreme assumptions such as mas- 
sive neutron stars of 2 Mq or more, which makes this 
not a very likely scenario. Furthermore, it seems ques- 
tionable whether the high entropies required to repro- 
duce the solar system r-process signature (e.g. Meyer & 
Brown 119971 Freiburghaus et al. I1999a|l can be provided 
by Type II supernovae (Witti et al. 119941 Qian & Woosley 
Hl&rJl Thompson et al. lSUUT)) . 

It should be noted that the r-process scenarios as- 
sociated with core-collapse supernovae discussed above 
are affected by considerable theoretical uncertainties since 
the physics of these events is not really understood (e.g. 
Rampp & Janka l2"UUm Liebendorfer et al. 200 1). Recently, 
Freiburghaus et al. l|1999bH presented for the first time 
detailed r-process calculations for neutron star mergers 
(NSM). Coalescing neutron stars potentially can provide 
in a natural way the large neutron fluxes required for the 
build-up of heavy elements through rapid neutron cap- 



ture. This scenario was subsequently explored and re- 
fined by Rosswog et al. ( 19991 l2"000(l . Their calculations 
showed that a few times 10~ 3 — 10 -2 Mq of r-process 
matter might be ejected in a merger event. This amount 
is significantly larger than the typical 10~ 5 — 10 -6 Mq 
of r-process material thought to be ejected in each core- 
collapse SN event (Cowan et al. HMTl Woosley et al. 11391 
Qian ESUTiOl Wanajo et al. l2lM|l . Since the rate of NSM in 
the Galaxy is significantly lower than that of Type II SNe 
(e.g. Tamman et al. 119941 Belczynski et al. I2002j) . either 
of these two sources may account for the total amount of 
r-process matter in the Galaxy (Qian & Woosley 119961 
Rosswog et al. 1199 9 ; Rosswog & Davies l2"002l Thielemann 

et ai.Enni. 

However, Qian (2000) argued that if NSM were a ma- 
jor r-process source, the low coalescence rate of binary 
compact objects would prevent any correlation between 
abundances of r-process elements and iron, which is clearly 
in disagreement with observed r-process abundances in 
stars more metal-rich than [Fe/H] ~ —2.5. In addition, 
the large amount of r-process ejecta from NSM would 
lead to r-process abundances in metal-poor halo stars that 
would be too large compared to observations. These ar- 
guments deserve closer examination, especially since no 
detailed chemical evolution calculations exist to date that 
assume NSM to be the major r-process source. Chemical 
evolution studies that treat SN with either lower-mass 
(8 — 10 Mq) or higher-mass (> 20 Mq) progenitors as the 
major r-process source were carried out by Ishimaru & 
Wanajo l|1999|) . Travaglio et al. I|1999|l and Tsujimoto et 
al. (2000). Unfortunately, each of these studies favoured 
different r-process sites: lower-mass SNe (Travaglio et al. 
1999), higher -mass SNe (Tsujimoto et al. U000), or both 
(Ishimaru & Wanai o |T9i39| are acceptable from the view- 
point of chemical evolution. 

The purpose of this work is twofold. First, neutron star 
mergers are included as major r-process site in detailed 
chemical evolution calculations for the first time. Second, 
we compare the r-process enrichment of the ISM under the 
assumption that neutron star mergers, lower-mass SNe II 
or higher-mass SNe II are the dominant r-process sites in 
the framework of inhomogeneous chemical evolution. In 
Sect. |21 the inhomogeneous chemical evolution model is 
presented. The impact of SNe II and NSM as r-process 
sites on Galactic chemical evolution is studied in Sect. |3] 
and conclusions are presented in Sect.QJ 

2. The chemical evolution model 

The large scatter seen in the abundances of neutron cap- 
ture elements in metal-poor halo stars is generally at- 
tributed to local chemical inhomogencities of the inter- 
stellar medium during the earliest stages of Galaxy for- 
mation and enrichment (e.g. McWilliam et al. (1995a 
I1995bl Ryan et al. HM® Burris et al. I2UUU1 Mishenina 
& Kovtvukh I2001fl . Several authors developed chemical 
evolution models (Ishimaru & Wanajo 119991 Raiteri et 
al. ITM91 Tsujimoto et al. ISUUUl Argast et al. I2UUU1 
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Travaglio et al. I2UUT1 Fields et al. l2T3U2jl that investigate 
the impact of local chemical inhomogeneities on the en- 
richment of the ISM and the scatter in element abun- 
dances. Although the methods employed by these authors 
differ significantly, the model results are in qualitative 
agreement. The notion of inhomogeneous chemical evo- 
lution therefore seems to be well established. 

In the present work we use the stochastic chemical evo- 
lution model presented in Argast et al. I|2000l hereafter 
Paper I) to investigate the enrichment of the ISM with 
r-process elements. In the following, we give a short sum- 
mary of the model and discuss some major changes (mass 
infall, star formation rate, SN la events and r-process 
sites) that were implemented for this work. 

2.1. Basic ingredients 

We model the chemical enrichment of the halo ISM within 
a volume of (2.0 kpc) 3 , down to a resolution of (50 pc) 3 . 
Primordial matter is assumed to be falling into this vol- 
ume, following an infall law of the form 



M(t) = a-t b -exp(-t/r). 



(1) 



Instead of specifying parameters a and b directly, we 
use the fact that with this description the time of maximal 
infall, t max , is given by 

tmax = b-T, (2) 

and that the total mass, M to t, falling into the volume is 



Mtat = 



a ■ t b ■ exp(— tjr) dt, 



(3) 



where t en d is the age of the system. The infall 
law is therefore fully described by the parameter set 

{Mtot: tmax: tend}- 

A crucial ingredient of chemical evolution models is 
the treatment of star formation (SF). Unfortunately, it is 
also one of the least understood. In this work we adopt 
a simple SF description based on Schmidt's law (Schmidt 
1959). The number n (t) of stars that are formed per time- 
step (of duration 10 6 yr) is determined by the product 



N t , 



n(t) 



N,, 



(4) 



i=l 



where pi (t) is the local ISM density at time t in a cell 
of volume (50 pc) 3 and the sum goes over all cells in the 
volume (N to t). The parameter v determines the star for- 
mation efficiency, and possible values for a range from 1 
(SF proportional to gas density) to 2 (SF triggered by 
cloud-cloud collisions, Larson H991f) . The factor N^ D \ was 
introduced to keep the number of newly formed stars in- 
dependent of the model resolution, so that Eq. 0] gives 
n(t) = v (p (t)) a if the ISM is homogeneously distributed 
(i.e. Pl (t) = (p(t)) for all*). 

The mass of a newly formed star is chosen randomly, 
subject to the condition that the mass distribution of all 



stars follows a Salpeter initial mass function (dN/dm oc 
rn~ 2 ' 35 ) with lower and upper mass limits of to; = 0.1 Mq 
and m up — 50 Mq , respectively. The average mass (to) of 
a star is given by the initial mass function (IMF) as 



(to) 



f mup to • dN 

J in 1 

f m " P dN ' 



(5) 



The number n (t) of star formation events translates into 
an average star formation rate (SFR) (ip (t)) at each time- 
step by multiplying the average stellar mass (to), i.e. 



(V> (t)) = n(t)-(m) 



Nt, 



N tot r m 

i=l J mi 



■ dN 



rm vp 



(6) 



The cells that undergo star formation are also cho- 
sen randomly, though the probability for a cell to get se- 
lected scales with its density. In this prescription, patches 
of denser material, e.g. in the neighbourhood of a super- 
nova remnant, are predominantly chosen for star forma- 
tion events. Note that a cell is allowed to form stars only 
if it contains at least 50 Mq of gas. This restriction is im- 
posed, so that SF is not biased towards low mass stars. 
With the imposed limit, a star formed in a cell can be 
of any mass in the range 0.1 — 50 Mq. Consequently, no 
SF will occur until enough material has fallen into the 
volume to exceed this limit. Due to the stochastic nature 
of our SF prescription, small deviations from the aver- 
age SFR in Eq. 0have to be expected at each time-step. 
Newborn stars inherit the abundance pattern of the ISM 
out of which they formed, thus carrying information about 
the chemical composition of the ISM at the place and time 
of their birth. 

Low mass stars (< 1 Mq) do not evolve significantly 
during the considered time but serve to lock up part 
of the total mass, thus affecting the abundances of el- 
ements with respect to hydrogen. Stars of intermediate 
mass (~ 1 — 10 Mq) return most of their mass at the 
end of their stellar lifetime, leaving a white dwarf as stel- 
lar remnant. Stars in the mass range of 10 — 50 Mq are 
assumed to explode as SNe II, polluting the neighbour- 
ing ISM with their highly enriched ejecta. Stellar life- 
times are taken from the Geneva Stellar Evolution and 
Nucleosynthesis Group (cf. Schaller et al. 119921 Schaereret 
al. H993al Schaerer et al. H993bl Charbonnel et al. HU9ljl . 
Stellar yields of O, Mg, Si, Ca and Fe are taken from 
Thielemann et al. H199ri[> and Nomoto et al. I|1997|> and 
are scaled according to Samland l|1998(l to account for the 
global chemical enrichment of the Galaxy. In particular, 
Fe yields are reduced by a factor of two. 

SN events pollute the neighbouring ISM with their 
nucleosynthesis products and sweep up the material in 
a spherical, chemically well mixed shell. Here, it is as- 
sumed that each SN pollutes « 5 x 10 4 M of ISM (Ryan 
et al. UnSSI Shigeyama & Tsuiimoto H^gjl . Stars which 
form out of material enriched by a single SN II show an 
element abundance pattern that is characteristic of the 
yields for the particular progenitor of this SN II. This 
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Table 1. Parameter values of the standard model. 



Parameter 


Description 


Value 


V 


modelled volume 


8 kpc 


Mtot 


total system mass 


10 8 Mq 


T 


infall decline timescale 


5 Gyr 


tmax 


time of maximal infall 


2 Gyr 


tend 


age of the system 


14 Gyr 


V 


SF efficiency 


15 


a 


exponent of SF law 


1.5 


Mio 


lower IMF mass limit 


0.1 Mq 


M up 


upper IMF mass limit 


50 Mq 


-PsNIa 


SN la probability 


6- 10" 3 



will lead to a large scatter in element abundance ratios 
([el/Fe]) as long as local inhomogeneities caused by SN II 
events dominate the halo ISM. As time progresses, super- 
nova remnants overlap and the abundance pattern in each 
cell approaches the average defined by SN II yield patterns 
for different progenitor masses and the IMF. This leads to 
a decrease of the scatter in element abundance ratios at 
later times. 

To determine which intermediate mass stars form 
Type la SN events, we adopt the following simple proce- 
dure: With probability PsNIa a newly formed intermediate 
mass star has a companion in the same mass range and one 
of these two stars will end its life as SN la. The mass of the 
companion (again in the range 1 — 10 Mq and following 
a Salpeter IMF) is determined randomly and the SN la 
event occurs after the less massive of the two stars enters 
the red giant phase. Although this procedure admittedly 
is rather simple, it has the advantage that the SN la fre- 
quency is determined by only one free parameter, namely 
the probability PsNIa- This parameter is chosen in such a 
way that the slope observed in [a/Fe] abundance ratios 
at [Fe/H] > — 1 is reproduced. We chose 10 Mq as upper 
mass limit for SN la progenitors to be consistent with the 
lower mass limit of SN II. Since we use such a simplified 
description to determine the occurrence of SN la, the up- 
per mass limit of SN la progenitors is not very relevant: 
Assuming a Salpeter IMF, more than 99% of the binary 
systems will consist of stars with masses < 8 Mq each. 
Finally, the yields of Type la SNe are taken from Iwamoto 
et al. lfT§fj51 Model CDD2). Note, that if revised, i.e. re- 
duced electron capture rates are applied (see Brachwitz et 
al. 2000), the yields of model CDD2 correspond to those of 
model WDD2, which are more accurate in terms of central 
conditions. 

In Table^ w e list the parameter values adopted for our 
standard model. 

2.2. Treatment of r-process sources in the model 

An ideal element to trace the r-process enrichment of the 
ISM is the pure r-process element Europium; approxi- 
mately 97% of the solar Eu abundance was produced in 
r-process events (Burris et al. I2000JI . Unfortunately, only 



a small sample of Eu abundances at very low metallicities 
(i.e. [Fe/H] < —2.5) are available to date. In order to trace 
the r-process enrichment at lower metallicities, the well 
studied element Barium is also used in this investigation. 
Ba abundances in stars are dominated by the s- process, 
and only w 15% of the solar Ba abundance was formed by 
rapid neutron capture (Burris et al. I2000f) . However, the 
r-process fraction [Ba r /Fe] of Ba abundances in stars can 
be easily calculated by subtracting the s-process compo- 
nent from its total Ba abundance (c.f. Burris et al. 2000). 
Since we did not include s-process sources in our model, 
displayed Ba r abundances in Sect. [21 of model and halo 
stars only show the r-process contribution to the total Ba 
abundances. 

In the following, r-process yields of Eu and Ba are es- 
timated under the assumption of a robust r-process for 
nuclei more massive than Ba, i.e. only one source is re- 
sponsible for the enrichment of the ISM with r-process el- 
ements beyond Z=56. In this work we focus only on these 
heavier neutron capture elements, since the enrichment of 
lighter neutron capture elements may require a second, 
independent source (e.g. Wasserburg et al. .ljJ9|J|. Thus, 
r-process matter ejected in an r-process event is consid- 
ered to consist of pure r-process nuclei more massive than 
Ba. First, the mass fractions of Ba and Eu are calculated 
with the help of Table 5 in Burris et al. ( 2000) and the 
standard solar system element abundances of Anders & 
Grevesse itTUSfl: 



(7) 



ATg a - < m Ba >_ _ ? ? , n - 2 

/Eu - ^ ATr ~ — « 1-2 ' 10 , 

22i N l- <mi> 



(8) 



where N[ is the number fraction of r-process nuclei of a 
nuclear species i, < m, > its mean atomic weight and the 
sum goes over all elements beyond Ba (Z > 56). Final 
yields are estimated by multiplying these mass fractions 
with the mass M r ejected in an r-process event, i.e. 



T Ba « 7.7 • 10~ 2 ■ M r , 
Feu ~ 1-2 • 10~ 2 ■ M r . 



(9) 
(10) 



The total ejected r-process matter, M r , is either deduced 
from observational evidence and/or theoretical consider- 
ations or is treated as a free parameter that is adjusted 
so that the results of the chemical evolution model are 
consistent with observations. 

We now shortly discuss the implementation of lower- 
mass SNe (8 - 10 Mq), higher-mass SNe (> 20 M ) and 
NSM as r-process sites in our model. We did not consider a 
hybrid model, where different r-process sites contribute to 
the enrichment of the ISM with neutron capture elements. 
Instead, it is assumed that for each of the three cases the 
corresponding site is responsible for the whole r-process 
element inventory (with Z > 56) of the Galaxy, i.e. is the 
major source of r-process elements beyond Ba. 



D. Argast et al.: Neutron star mergers vs. SNe II as dominant r-process sites 



5 



2.2.1. r-process yields from core-collapse SNe in the 
mass range 8 — 10 Mq 

Lower mass core-collapse SNe (8 — 10 Mq) are suspected 
to be major r-process sources without contributing much 
iron-peak elements to the enrichment of the ISM (Wheeler 
et al. H998)l . Recently, Sumiyoshi et al. |2001| presented re- 
sults of r-process nucleosynthesis occurring during the ar- 
tificially induced prompt explosion of an 11 Mq star. The 
resulting distribution of r-process yields is in reasonable 
agreement with solar r-process abundances for the heav- 
ier elements, whereas lighter r-process elements (A < 100) 
are underproduced. Since elements in this mass-regime are 
often overproduced in neutrino driven wind models, the 
authors speculate that lower-mass SNe II might be re- 
sponsible for the production of massive r-process nuclei, 
whereas higher-mass SNe II account for the less massive r- 
process nuclei. Note however, that it is still doubtful if the 
prompt explosion of a massive star may occur in reality 
(Liebendorfer et al. 2001; Sumiyoshi et al. 2001). Chemical 
evolution models including lower mass SNe as r-process 
source were presented by Ishimaru & Wanajo fl222$ and 
Travaglio et al. l(T^jl . 

In the treatment of such lower-mass SNe, we proceed 
similar to Ishimaru & Wanajo 1)1999(1 : In model SN810, 
r-process nucleosynthesis is assumed to occur in the mass 
range 8 — 10 Mq with constant Ba and Eu yields over 
the whole mass range and assuming a Salpeter IMF. The 
yields then are deduced from the average [Eu/Fe] and 
[Ba'/Fe] ratios (both i=s 0.5) of metal-poor halo stars: 
F B a = 5.3-10- 7 M Q andy Eu = 8.3-10" 8 M Q , which equals 
to a total of w 6.9 • 10~ 6 Mq of ejected r-process matter 
per event (Eas. I§1 and 110)1 . Furthermore, we assume that 
the amount of a- and iron peak elements synthesized in 
these SN II events are negligible (Hillebrandt et al. 11984)1 . 

2.2.2. r-process yields from core-collapse SNe more 
massive than 20 Mq 

Nucleosynthesis of r-process elements might also occur 
in neutrino driven winds or jets from nascent neutron- 
stars during the delayed explosion of high-mass stars 
with masses > 20 Mq (e.g. Woosley k Hoffman ITM?! 
Takahashi et al. IT^Ml Woosley et al. HUMl Qian & 
Woosley H^HI Thompson et al. I2UUT1 Nagataki IMiTl 
Wanajo et al. EMI EES Terasawaet al. lMT^I . The phys- 
ical conditions needed for a robust r-process in neutrino 
driven winds (such as high entropies, low electron fraction 
and/or short dynamical timescales) are hard to achieve 
in present models, and lighter r-process nuclei are often 
overproduced (e.g. Thompson et al. 12001)1 . As an exam- 
ple, Wanajo et al. ( 2002) require a very massive and com- 
pact proto-neutron star of 2.0 Mq to reproduce the so- 
lar system r-process abundance pattern of heavier neu- 
tron capture elements. By allowing for strong magnetic 
fields, Thompson (2003) finds that the required physical 
conditions are much more easily achieved than in previ- 
ous models. Yet, there exist still some considerable theo- 



retical uncertainties and only very high-mass stars seem 
to be capable of providing the environment for a ro- 
bust r-process which reproduces the abundance pattern 
of neutron-capture elements beyond Ba. 

Williams (1987) reported the detection of Ba absorp- 
tion lines in the spectra of SN 1987A (progenitor mass 
f» 20 Mq), and Mazzali et al. 1)1992)1 noted that Ba 
was lacking at the very surface of the ejecta. Because of 
this remarkable feature, Tsujimoto & Shigeyama (2001) 
identify SN 1987A as r-process site and estimate that 
w 6 ■ 10~ 6 Mq of Ba were synthesized during the SN 
event. However, the poor knowledge of the far UV ra- 
diation field in the envelope of SN 1987A makes it dif- 
ficult to determine a truly reliable Ba abundance in its 
ejecta (Utrobin & Chugai 2o03 and it is not yet estab- 
lished without doubt that the detected Ba was formed in 
the SN event. Complementary to the observation of Ba 
in SN 1987A, Tsujimoto et al. ()2"UUU|> deduce Ba and Eu 
yields from observations of metal-poor halo stars and in- 
homogeneous chemical evolution models, suggesting that 
core-collapse SNe in the mass range 20 — 25 Mq dominate 
the production of r-process elements. They propose a Ba 
yield of 8.5 • 10~ 6 Mq for a 20 Mq and 4.5 • 10~ 8 M Q for 
a 25 Mq progenitor and a Eu yield of 1.3 • 10 _6 Mq and 
7.0 • 10 -9 Mq for a 20 and 25 Mq progenitor, respectively. 

In a first model (SN2025), we adopt yields similar to 
the ones given by Tsujimoto et al. 1)2000)1 and assume that 
SNe outside the indicated mass range (i.e. 20 — 25 Mq) do 
not contribute significantly to r-process nucleosynthesis. 
In this model, a 20 Mq star produces 4.3 • 10~ 6 Mq of 
Ba and 6.5 ■ 10 -7 Mq of Eu, whereas a 25 Mq star yields 
2.3 • 10" 8 Mq of Ba and 3.5 • 10~ 9 M Q of Eu. For stars 
between 20 and 25 Mq, the yields are interpolated linearly. 
Thus, a SN II in this mass range ejects on average « 3 • 
10 -5 Mq of r-process matter. Compared to the r-process 
yields proposed by Tsujimoto et al. 1)2000)1 our yields are 
reduced by a factor of two. This is the case, since we had 
to scale the Fe yields of Thielemann et al. ( 199B|) by the 
same factor to account for the global chemical enrichment 
of the Galaxy. 

As an alternative to model SN2025, we also calculated 
the chemical evolution of the ISM with r-process yields 
from SNe II in the mass range 20-50 Mq (model SN2050). 
The yields were chosen in such a way that the whole range 
of r-process abundances in metal-poor halo stars are re- 
produced in the model. The r-process yield of a 20 Mq 
star was set to 1.5 • 10~ 4 , which is about three times the 
r-process yield of a 20 Mq star in model SN2025. This 
large amount of ejected r-process matter requires a rapid 
decline in the r-process yields of progenitors with masses 
20 Mq < to < 28 Mq. Otherwise, the average [r/Fe] ra- 
tios of model stars do not reproduce the observed average 
abundances of metal-poor halo stars ([r/Fe] « 0.4 — 0.5). 
Yields of progenitors in the mass range 28 — 50 Mq then 
are assumed to be constant. On average, 1.2 • 10~ 5 Mq of 
r-process matter are ejected in each event. 

Table El lists our Ba yields Yea (m) and Eu yields 
Yeu (m) and total ejected r-process matter M r as func- 



6 



D. Argast et al.: Neutron star mergers vs. SNe II as dominant r-process sites 



Table 2. Adopted Ba (Y Ba (to)) and Eu (Y Eu (m)) yields 
and ejected r-process matter (M r ) as function of progeni- 
tor mass m of models SN810, SN2025 and SN2050. Yields 
in the mass range 8 — 10 Mq are assumed to be constant 
(SN810). Yields in the mass range 20 — 25 Mq are linearly 
interpolated (SN2025). Yields of model SN2050 exhibit 
a more complex behaviour: A rapid decline in the mass 
range 20 — 28 Mq is assumed. Progenitors with masses 
28 — 50 Mq all have the same, constant yield. 



Table 3. Constraints on the Galactic NSM rate /nsm and 
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Remarks: "Model SN810, "Model SN2025, c Model SN2050 

tion of progenitor mass m of models SN810, SN2025 and 
SN2050. We point out that the yields adopted in this work 
are chosen in such a way that observational constraints 
from metal-poor halo stars are respected. Whether these 
yields are feasible would have to be tested with the help 
of self-consistent core-collapse SN models, which, unfortu- 
nately, are not available to date. 

2.2.3. r-process yields from neutron star mergers 

Another major source of r-process elements might be neu- 
tron star mergers (Freiburghaus et al. I1999bl Rosswog et 
al.^999 2000]). A scheme similar to the determination of 
SN la events was applied for the determination of their oc- 
currence (c.f. Sect. 12. ljl : With probability Pnsm a newly 
formed high mass star (in the range 10 — 50 Mq) has a 
companion in the same mass range. It is assumed that, 
subsequent to both SN II events, the remaining neutron 
stars will eventually merge. The time t c that is needed for 
the coalescence of the two neutron stars and the proba- 
bility Pnsm for the occurrence of NSM events are treated 
as free parameters which can be adjusted within given 
observational estimates on the coalescence rate of binary 
compact objects and merger timescales (van den Heuvel 
& Lorimcr 1996. Kalogera & Lorimer 2000; Belczynski et 
al. 12)02). As was the case for SN la events, this treatment 
is very simple, but has the advantage that all difficulties 
associated with the determination of formation channels 
of neutron star binaries (or neutron star / black hole bina- 
ries) and the corresponding formation rates are combined 
in only one free parameter. Furthermore, the lower mass 



ejected r-process matter M r . 
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boundary for neutron star formation (here 10 Mq) is of 
little importance: The decisive term for the enrichment 
of the ISM with r-process material in our model is the 
NSM rate. Since a decrease of the lower mass boundary 
results in an increase of the NSM rate, a simple adjust- 
ment of the parameter Pnsm is sufficient to retrieve a NSM 
rate consistent with constraints on the galactic rate. In 
the model the ejected r-process matter is also mixed with 
en 5 • 10 4 Mq of ISM and is assumed to be distributed in a 
spherical, chemically well mixed shell. This might not be 
true, however, since it is well conceivable that r-process 
matter is ejected only in the orbital plane of the merging 
neutron stars. In this case, the ejecta would be distributed 
over a smaller volume and would consequently mix with 
a smaller amount of ISM than the one assumed in our 
model. As will be seen below (Sect. ET2l . this would lead 
to even higher [r/Fe] ratios in the computed model stars. 
Thus, our assumptions are in favour of the NSM scenario 
and changing the mixing mechanism in our model would 
only aggravate the problems associated with the high r- 
process yields from NSM. 

Coalescence timescales for neutron star mergers are 
typically estimated to be of the order 100 — 1000 Myr (e.g. 
Portegies Zwart & Yungel'son ITMSI Fryer et al. ESS). 
Recently, Belczynski et al. i|200 2) suggested a dominat- 
ing population of short lived neutron star binaries with 
merger times less than 1 Myr. This population of neutron 
star mergers might be formed through channels involv- 
ing mass-transfer episodes from helium stars, leading to 
tightly bound binary systems with very short orbit de- 
cay timescales. However, this result depends on the oc- 
currence of a common envelope (CE) phase of the progen- 
itor He stars, which was treated in a simplified manner by 
Belczynski et al. f2002)| . The authors caution that detailed 
hydrodynamical calculations of the CE phase still have to 
confirm these results. In view of these uncertainties, we 
adopt two different coalescence timescales t Cl namely 1 and 
100 Myr. In each case, it is assumed that all neutron star 
binaries merge within this time. Coalescence timescales 
of the order 1 Gyr have not been considered since they 
are not consistent with observations of neutron capture 
elements in ultra metal-poor stars (c.f. Sect EH|) . 

Constraints on the Galactic NSM rate are also con- 
troversial. Van den Heuvel & Lorimer (1996) estimate a 
rate of approximately 8 • 10~ 6 yr -1 , Kalogera & Lorimer 
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Table 4. Parameter values adopted for NSM events (see 
text for details). For each value of Pnsm, two models with 
the indicated coalescence timescales t c have been calcu- 
lated. 
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( 2000) give an upper limit of (7.5 • 10~ 7 - 1.5 • 10" 5 ) yr" 1 
whereas Belczynski et al. I|2002[) get rates in the range 
(l0~ 6 — 3 • 10~ 4 ) yr . In Newtonian calculations, the 
amount M r of r-process matter ejected in a NSM event is 
of the order of a few times 10~ 3 — 10~ 2 Mq, depending on 
the initial configuration of the binary system (Rosswog et 
al. I1999f) . Taking general relativistic effects into account, 
Oechslin et al. ljM%)> find 5 • 10~ 5 - 2 • 10~ 4 M as lower 
limit for M r . All constraints are again listed in Table |3| 

Table0]lists the parameter values adopted for the NSM 
probability Pnsm, the resulting average NSM rate /nsm, 
the amount of r-process matter M r ejected in each event 
and coalescence timescales t c . After a value for Pnsm has 
been chosen, the coalescence timescale is the only free pa- 
rameter remaining in the model since the actual NSM rate 
is determined by the NSM probability and the SF rate (re- 
sulting from the parameter values given in Tabled . In ad- 
dition, the NSM rate and the amount of ejected r-process 
matter are tightly correlated, since the total amount of 
r-process matter in the Galaxy (w 10 4 Mq, Wallerstein et 
al. I1997f) has to be reproduced. Thus, higher NSM rates 
require that less r-process matter is ejected in each event, 
and vice versa. 

Comparing the values in Table 01 with the constraints 
listed in Table El reveals that the NSM rate of the first row 
in Table 0] is too high by a factor of 10-100 and that the 
corresponding low value of 10~ 4 Mq of ejected r-process 
matter is close to the lower limit given by Oechslin et al. 
(2002). NSM rates listed in the three lower rows seem to 
be consistent with the galactic NSM rate. However, the 
ejected r-process matter in the last row is clearly beyond 
the upper limit of allowed values. 

3. Enrichment of the ISM with r-process elements 

3.1. SN II as dominating r-process sites 

In this section, the enrichment of the ISM with neu- 
tron capture elements is discussed under the assump- 
tion, that the dominating r-process sources are either 
lower-mass SNe II (8 — 10 Mq) or higher-mass SNe II 
(> 20 M Q ). The results of models SN810, SN2025 and 
SN2050 are shown in Figs.n~El respectively. The figures 
show the evolution of [Eu/Fe] and [Ba r /Fe] as function of 
metallicity [Fe/H]. Model stars are shown as black dots, 
whereas observations are indicated by filled red squares 
and diamonds. Observations are taken from Peterson et 



al. <HHn|), Gratton & Sneden p991al IT99TB|) . Ryan et al. 
(HniIJ), Edvardsson et al. l(T^9ljl . Francois et al. iTTSgljl . 
Beveridge & Sneden (TH?9l|l . McWilliam et al. I|1995afl . 
Ryan et al. (JEHU), Jehin et al. lfTM5| . Aoki et al. |2T)0T)| . 
Burris et al. |2000 ), Mashonkina & Gehren p)00f> . Sneden 
et al. I|2000a|l . Mashonkina & Gehren l|2001|l . Mishenina 
& Kovtyukh (HHU, Koch & Edvardsson pUU2|l and 
Stephens & Boesgaard ( 2002 ) . Single observations of stars 
are marked by a square. In the case where multiple obser- 
vations of a star are present, we plot the most recent one if 
all observations were published before the year 2000 (also 
marked by squares). If several more recent observations 
are available, the given element abundances are averaged 
(marked by diamonds). The yellow line shows the average 
element abundances in the model ISM and can directly be 
compared to classical chemical evolution models, which 
assume that the ISM is well mixed at all times. Purple 
filled circles with error bars denote the average [r/Fe] ra- 
tios of model stars and their standard deviation in [Fe/H] 
bins with binsize 0.1 dex. Note, that the [Ba r /Fe] plots 
only show the r-process contribution to the total Ba abun- 
dances of halo stars. According to Burris et al. I|2000[) . r- 
process Ba abundances can be computed by removing the 
s-process contribution to Ba in stars with [Fc/H] > —2.5, 
if Eu abundances have also been determined. For stars 
with metallicities [Fe/H] < —2.5, it can be assumed that 
the whole Ba inventory is of pure r-process origin. Thus, 
published Ba abundances of such metal-poor halo stars 
simply have to be renormalized to the level of the solar 
r-process fraction (i.e. «s 15% of the total solar Ba abun- 
dance). Care has been taken to remove known carbon stars 
from our sample. Such stars mostly show unusually large 
Ba abundances, which arc thought to originate from mass 
transfer of s-process enriched matter in binary systems. 

The evolution of r-process elements shown in Figs.^ 
13 are qualitatively very similar. At very low metallicities 
([Fe/H] < —2.5), a large scatter in abundances of model 
stars is visible. This scatter is due to chemical inhomo- 
geneities in the early ISM (c.f. Paper I). The scatter de- 
creases as the mixing of the ISM improves and finally 
reaches the IMF averaged mean. At this stage, the ISM 
can be considered well mixed and the further evolution is 
comparable to the one of classic chemical evolution mod- 
els. In the following, we point out some important fea- 
tures of the ISM enrichment resulting from models SN810, 
SN2025 and SN2050: 

1. In all models, r-process elements appear very early in 
the enrichment of the ISM. Some model stars with r- 
process abundances exist even at metallicities [Fe/H] 
< —4, which is in agreement with observed Ba abun- 
dances in ultra metal-poor stars. Neutron capture el- 
ements at such low metallicities are considered to be 
of pure r-process origin, since s-process enrichment is 
expected to occur at a later stage in the chemical evo- 
lution of the Galaxy (e.g. Travaglio et al. I1999|l . 

2. For metallicities [Fe/H] > —2, the scatter in [Eu/Fe] 
and [Ba r /Fe] abundances of model stars is compa- 
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Fig. 1. Evolution of [Eu/Fe] and [Ba r /Fe] abundances as function of metallicity [Fe/H]. Lower-mass SNe II (8 — 10 Mq, 
Model SN810) are assumed to be the dominating r-process sources. Black dots denote model stars, observations are 
marked by filled squares and diamonds (see text). Average ISM abundances are marked by a continuous line. Filled 
circles with error bars denote average abundances of model stars and their standard deviation in [Fe/H] bins with 
binsize 0.1 dex. 




Fig. 2. Evolution of [Eu/Fe] and [Ba r /Fe] abundances as function of metallicity [Fe/H]. Higher-mass SNe II (20 — 
25 Mq, Model SN2025) are assumed to be the dominating r-process sources. Symbols are as in Fig. ^ 



rable to the one seen in observations. The main dif- 
ference between the higher-mass (models SN2025 and 
SN2050, Fig. H and and lower-mass (model SN810, 
Fig-CJ SNe II scenarios is the larger scatter in [r/Fe] at 
[Fc/H] « — 3 and the slightly more pronounced scatter 
at higher metallicities of models SN2025 and SN2050. 
One might argue, that the model stars in Figs. [21 and 



El at [Fe/H] > -2 with high [r/Fe] abundances (> 1) 
should also be visible in the observational data. The 
bulk of model stars, however, has abundances [r/Fe] 
w 0.4 — 0.5 and it would thus be very unlikely to ob- 
serve stars with high [r/Fe] ratios in this metallicity 
regime. The situation below [Fe/H] = —2 is differ- 
ent: Only few model stars are present in this metallic- 
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Fig. 3. Evolution of [Eu/Fe] and [Ba r /Fe] abundances as function of metallicity [Fe/H]. Higher-mass SNe II (20 — 
50 Mg, Model SN2050) are assumed to be the dominating r-process sources. Symbols are as in Fig. ^ 



ity range and their [r/Fe] abundances are evenly dis- 
tributed and do not clump around the IMF averaged 
value. 

3. Models SN810 and SN2025 fail to reproduce the two 
metal-poor stars with the highest Eu abundances 
([Eu/Fe] > +2.0 at [Fe/H] w -2.5). The existence 
of such ultra r-process enhanced stars might pose a 
serious problem for the SN scenario, since an unusual 
large amount of ejected r-process matter (> 10 _4 Mq) 
is required to reproduce these observations in chemi- 
cal evolution models (model SN2050, Fig. EJ). In the 
neutrino driven wind model of Wanajo et al. (2001), 
massive proto-neutron stars with a high compaction 
ratio are required to achieve r-process yields of this or- 
der of magnitude. However, chemical evolution models 
constrain such large r-process yields to only a small 
progenitor mass range (see Table [5} ■ Otherwise, the 
total r-process inventory of the Galaxy is exceeded. 
On the other hand, Qian & Wasserburg PUUTI ISUOljl 
argue that the highest Eu and Ba abundances in metal- 
poor halo stars may originate from the dumping of r- 
process enriched matter from a high-mass companion 
onto the surface of a low mass star. This explanation 
requires that during the SN II event of the high-mass 
companion, r-process elements but no iron-peak nuclei 
are produced. Possible candidates for this kind of SN II 
events are progenitors in the mass range 8 — 10 Mq; 
see e.g. Qian & Wasserburg (2002) for a discussion of 
this mechanism. To date, our model is not able to han- 
dle the occurrence of surface contamination in binary 
systems, so that the results of model SN810 do not 
directly apply to this scenario. We speculate however, 
that if this mechanism would be at work during the 



evolution of the Galaxy, we might also expect stars 
with high [r/Fc] abundances at metallicities [Fe/H] 
> —2. No such stars have been observed to date. Yet, 
since many abundance studies primarily target metal- 
poor stars, it is possible that r-process enriched, metal- 
rich stars simply have been missed. Note, that if iron 
is produced in the SN II event of the companion, the 
resulting surface abundance will be a mixture of the 
intrinsic r-process abundances and the yields of the 
SN II. In this case, large r-process yields are still re- 
quired to enrich the companion to the high level seen 
in observations. 

4. The models also predicts some stars with very low 
[Eu/Fe] and [Ba r /Fe] ratios at [Fe/H] < -2, which 
are not observed to date. However, if stars with [r/Fe] 
< — 1 exist in the galactic halo, their Eu and Ba lines 
may be too weak to be detectable (c.f. observational 
limits quoted in Travaglio et al. (J2001)). In the mod- 
els, stars with very low [Eu/Fe] ratios are inevitably 
produced if only a limited mass-range of SNe II con- 
tributes to the r-process enrichment. In this case, pock- 
ets without or with only small amounts of r-process 
elements may form in the halo ISM if primarily SNe II 
without r-process yields contributed to the local en- 
richment. 

5. Model stars in Fig. ^ with high r-process abundances 
([r/Fe] > +2) at metallicities below [Fe/H] < -3.5 
owe their origin to the same phenomenon. Here, lower- 
mass SNe II are the primary source of r-process ele- 
ments, without contributing much iron-peak nuclei to 
the enrichment of the ISM. Thus, pockets with large r- 
process abundances may form during the early enrich- 
ment of the ISM, similar to pockets with large Fe abun- 
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dances. However, [r/Fe] ratios soon converge to the av- 
erage [r/Fe] ratios of metal-poor halo stars and from 
then on give a good fit to observations. It is interesting 
to note, that this behaviour is typical for inhomoge- 
neous chemical evolution models and is not present in 
classic chemical evolution models which assume that 
the ISM is well mixed at all times (e.g. Travaglio et 
al. I1999f ). This simplifying assumption leads to a slow 
increase of the average ISM [r/Fe] abundances, since 
r-process yields from 8 — 10 Mq are instantaneously 
mixed with iron yields from 10 — 50 Mq. 
We briefly note here, that the fit of model SN810 to ob- 
servations of metal-poor halo stars can be improved, if 
infall of pre-enriched ISM with [Fe/H] — — 3 is assumed 
(e.g. Qian & Wasserburg 120021) . In this case, model 
stars with high [r/Fe] abundances at [Fe/H] ~ —4 are 
shifted to [Fe/H] = —3, where such r-process enriched 
stars are observed. This is the case only for model 
SN810. The fit to observations of models SN2025, 
SN2050 and the NSM models (see Sect. E21> are not 
changed significantly by the assumption of metal-rich 
infall. 

6. The evolution of average [r/Fe] abundances of model 
stars (purple circles) as function of mctallicity is an 
important tool to compare model results with observa- 
tions: Slopes in the distribution of model stars, which 
might be hidden in the full [r/Fe] plots, will show up 
in the distribution of these averaged abundances and 
may be compared to trends seen in observational data. 
Note, that the evolution of average ISM abundances 
(yellow line) is generally different from the evolution 
of average model star abundances. This is the case, 
since ISM abundances are averaged over all cells in our 
volume, which is equivalent to assuming a well mixed 
ISM at all times and thus may be compared to classic 
chemical evolution models. Average model star abun- 
dances, on the other hand, are still affected by chemical 
inhomogeneities, since they depend on the number of 
stars with a given [r/Fe] abundance. 
This is demonstrated in Fig. ^ (model SN810). 
Although the calculated standard deviation of model 
stars is large (due to the low number of model stars 
with metallicities below [Fe/H] < —3), the average 
model star abundances rise with decreasing metallic- 
ity Thus, model SN810 would predict a rising slope in 
[Ba' /Fe], contrary to the one present in observations 
of metal-poor halo stars. In model SN2025, the average 
model abundances stay more or less constant, whereas 
in model SN2050, average abundances give a nice fit 
to the observational data. Note, that the empty bin at 
[Fe/H] = —4 in Fig.|5]only contains one model star. 

In summary, core-collapse SNe seem to be a valid 
source of r-process elements from the point of view of 
chemical evolution, since the enrichment of the ISM in 
the cases discussed above is in qualitative, if not necessar- 
ily quantitative, agreement with observations. Based on 
the evolution of average model star abundances and the 



range of the scatter at [Fe/H] > — 3, model SN2050 gives 
the best fit to observations, whereas model SN810 gives 
the worst. Note however, that the initial scatter in [r/Fe] 
ratios strongly depends on the yields as function of pro- 
genitor mass. As was discussed in Argast et al. H2002J) . the 
progenitor mass dependence of stellar Fe yields are not 
known to date and the distribution of observed element 
abundances as function of metallicity can easily be repro- 
duced by the appropriate choice of Fe yields. Additionally, 
r-process yields in this work are chosen completely ad hoc 
and yet lack any theoretical backup. In view of these un- 
certainties, we conclude in accordance with Ishimaru & 
Wanajo (199SJ), that it is not possible to date to rule out 
either lower-mass or higher-mass SNe II within the frame- 
work of inhomogeneous chemical evolution. 

3.2. NSM as dominating r-process sites 

The enrichment of the ISM with r-process elements in the 
case of NSM acting as major r-process sources is discussed 
in the following. As already mentioned in Sect. l2.2!^l eight 
model-runs were carried out in total, assuming NSM rates 
ranging from 2 ■ 10 _3 yr _1 to 2 • 10~ 6 yr~ 2 and coales- 
cence time-scales of 1 and 100 Myr (see Table 0J. A con- 
densed overview of all models displaying the evolution of 
[Ba r /Fe] is shown in Figs. 0] and [5] The evolution of r- 
process abundances is strikingly different from the case 
in which r-process nucleosynthesis occurs in SNe II (see 
Figs |3Jl • All models with NSM as r-process sources fail 
to reproduce observations. 

As a representative example, the model shown in the 
upper right panel of Fig. 01 is discussed in the following. 
It was calculated with a NSM rate of 2 ■ 10~ 4 yr _1 , a 
total of 10 -3 Mq of ejected r-process matter and a co- 
alescence time-scale t c of 1 Myr. Note, that the case dis- 
cussed here assumes a NSM rate which is at the upper 
limit given by present estimates of the galactic NSM rate: 
Estimates range from 8 • 10 _6 yr~ 1 (van den Heuvel & 
Lorimer lT3M| to (l0~ 6 - 3 • 10" 4 ) yr" 1 (Belczynski et al. 
2002 ) and Kalogera & Lorimer ( 2000 ) even give an upper 
limit of (7.5 • 10~ 7 — 1.5 • 10~ 5 ) yr -1 . Furthermore, the 
coalescence timescale adopted in this model is only 1 Myr, 
which is two to three orders of magnitude lower than the 
classical estimate of 100 — 1000 Myr (e.g. Portegies Zwart 
& Yungel'son 19iSI Fryer et al.[l999). However, Belczynski 
et al. (12002) suggest the existence of a dominating pop- 
ulation of short lived neutron star binaries with merger 
times less than 1 Myr. 

The following qualitative differences are immediately 
visible, if the upper right panel of Fig. 0] is compared to 
Figs. El 

1. r-process nuclei appear at the earliest around [Fe/H] 
pa —2.5, whereas r-process elements such as Eu are 
observed at [Fc/H] = — 3 and probably even down to 
[Fe/H] = —4 in the case of Ba r . The reason for this late 
injection of r-process matter in the model is the low 
NSM rate (compared to the SN II rate). The decisive 
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Fig. 4. [Ba r /Fe] vs. [Fe/H] for the NSM rates 2 • lO^yr -1 , 2 • lO^yr -1 , 2 • lO^yr" 1 and 2 • lO^yr -1 (from left to 
right and top to bottom). The coalescence timescale adopted in these cases is t c = 10 6 yr. 



parameter to obtain a given NSM rate in the model 
is the probability Pnsm which, in this case, is set to 
4 ■ 10" 3 (see Table |]J. Thus, on average « 250 SNe II 
events occur before the first NSM event and r-process 
nuclei are injected into an ISM already pre-enriched by 
SNe II. The lower the NSM rate, i.e. the lower P N sm, 
the later (in time and metallicity) the occurrence of 
r-process elements. This late injection of r-process nu- 
clei by NSM is the reason why we did not consider 
coalescence timescales of the order of 1 Gyr: In our 
model, 1 Gyr after the first SF event the metallicity 
of the ISM reached [Fe/H] w -0.9. The first appear- 



ance of r-process nuclei at this metallicity is clearly 
not consistent with observations. Thus, the advantage 
that the NSM scenario can produce a large scatter in 
[r/Fe] close to 3 dex, as observed in ultra metal-poor 
stars, turns into a disadvantage because it appears at 
too high metallicities where observations do not show 
this effect anymore. 
2. There is a prominent tail of model stars with very 
low [r/Fe] ratios at [Fe/H] = —2. Such model stars 
are also present in Figs. ^ - |21 In this case, how- 
ever, a significant fraction of them shows this prop- 
erty. Furthermore, the tail develops not until [Fe/H] 
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[Fe/H] [Fe/H] 



Fig. 5. Same as Fig.0]but with t c = 10 s yr. 



> —2, whereas in the SN II case there are virtually 
no model stars with very low [r/Fe] ratios above this 
metallicity. Since the tail develops at higher metallic- 
ities, it can not be dismissed as unobservable. Even if 
r-process element abundances with [r/Fe] < — 1 could 
not be detected, the upper part of the tail should be 
visible in the observational data. 

3. Even at late times in the enrichment of the ISM ([Fe/H] 

> —1), the scatter in possible [r/Fe] ratios is of the 
order 1.5 — 2.0 dex, whereas observations of [Eu/Fe] 
and [Ba r /Fe] abundances show a scatter of approxi- 
mately 0.2 — 0.3 dex. The large scatter at high [Fe/H] 



occurs only for NSM but not for SNe II. Since the 
total amount of r-process matter ejected in a sin- 
gle NSM event (here 10 _3 Mq) is very large (cf. to 
rj 10~ 6 - 10~ 5 Mq for SNe II) and the frequency of 
NSM is very low (rj 100 - 1000 times lower than that 
of SNe II), NSM events may still cause significant lo- 
cal chemical inhomogeneities, in spite of the advanced 
enrichment of the ISM by many SNe II. 

These aspects strongly argue against NSM as the dom- 
inating r-process source, especially since the parameters 
used for the model in discussion are at the upper limit set 
by theoretical and observational constraints. Lower NSM 
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Fig. 6. loge (Ba r ) vs. [Fe/H] of model SN2050 (left panel) and NSM occurring at a rate of 2 • 1(T 5 yr" 1 (right panel). 
Symbols are as in Fig. 2] Additionally, numbered circles in the left panel give loge (Ba r ) of SNe II with the indicated 
progenitor mass, when mixed with 5 • 10 4 Mq of primordial ISM. Similarly, the labelled circle in the right panel 
indicates loge (Ba r ) of NSM when mixed with 5 • 1O 4 M of enriched ISM (at [Fe/H] = -2). 



rates and, consequently, larger ejecta masses, only aggra- 
vate the problems mentioned above. The dramatic changes 
in the distribution of r-process abundances occurring with 
decreasing NSM rate are clearly visible in the sequence of 
plots in Figs. 0] and The NSM rates adopted in the mod- 
els decrease from left to right and top to bottom by one 
order of magnitude for each panel, i.e. from 2 ■ 10 _3 yr _1 
for the upper left panel down to 2 • 10~ 6 yr _1 for the lower 
right panel. Figs. 0] and assume coalescence timescales 
of 1 and 100 Myr, respectively. 

In the upper left panels of both figures (NSM rate 
2 • 10~ 3 yr _1 ), a few model stars with r-process abun- 
dances first appear around [Fe/H] w —3.5, in contrast to 
[Fe/H] « —1 in the panels at the lower right (NSM rate 
2 • 10 _6 yr _1 ). Simultaneously, the scatter in [Ba r /Fe] and 
[Eu/Fe] at solar metallicity, which is of the same order 
as the one observed in the upper left panels, increases to 
almost 3 dex in the lower right panels, clearly not con- 
sistent with observations. Furthermore, the tail of model 
stars with low [r/Fe] abundances gets more pronounced 
and concentrated and is shifted to higher metallicities. 

Contrary to expectations, the impact of the coales- 
cence timescale t c on the distribution of r-process nuclei 
is rather small. This is surprising, since a large value of t c 
results in a delayed injection of r-process matter into the 
ISM. However, this can be understood by noting that in 
the NSM models, the enrichment of the ISM was adjusted 
to proceed very slowly at the beginning of Galaxy forma- 
tion: It takes w 200 Myr to enrich the ISM from [Fe/H] 
= —4 to [Fe/H] = —3. Thus, a coalescence timescale of 100 
Myr is of the same order as the halo enrichment timescale, 



resulting in only a slight shift of the first r-process mat- 
ter injection to higher metallicities. The slow enrichment 
was assumed in order to enable r-process nuclei injection 
from NSM even at low metallicities. However, as Figs. 0] 
and [3 show, this is not the case even with the help of 
a slow enrichment history. Thus, a fast halo enrichment 
would shift the moment of first r-process matter injection 
to even higher metallicities. Note that, on the basis of 
these considerations, coalescence timescales of the order 
1 Gyr clearly are not consistent with the occurrence of 
r-process nuclei early on during Galaxy formation. 

Fig. demonstrates another qualitative difference be- 
tween the SN II and NSM scenarios. In this figure, 
loge(Ba r ) 1 vs. [Fe/H] of model SN2050 (left panel) and 
NSM occurring at a rate of 2 • 10 -5 yr _1 (right panel) are 
shown. In the SN II scenario, the scatter in loge (Ba r ) at 
low metallicities ([Fe/H] < —2) is primarily determined by 
the progenitor mass dependent r-process yields of SNe II 
(indicated by numbered circles). Model stars at [Fe/H] 
< —2 with Ba r abundances above loge(Ba r ) = —3 are 
formed out of material enriched by a single SN II whereas 
some model stars with loge(Ba r ) < —3 are formed out 
of material which underwent subsequent mixing with pri- 
mordial ISM. Stars with such low r-process abundances 
disappear at [Fe/H] > —2 due to the increased enrichment 
of the ISM with r-process material. 

In the NSM case, the scatter in loge(Ba r ) is solely 
determined by subsequent mixing of r-process enriched 
material with r-process poor ISM. Under the assumption 
that a NSM event is able to enrich 5 • 10 4 M of ISM with 

1 loge (A) = log (7V A /iVff) + 12 
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10 -2 Mq of r-process matter, model stars with typically 
loge (Ba r ) = 2 or more can be formed (indicated by the 
labelled circle in the right-hand panel of Fig. |BJ). Thus, 
the distinct peaks in the distribution of Ba r abundances 
of model stars are telltale signs of recent NSM events. 
However, since NSM events are rare, there is time enough 
to thoroughly mix the r-process enriched material in the 
vicinity of a recent NSM with r-process poor ISM. In con- 
sequence of this mixing and the low NSM rate, the bulk 
of model stars with loge (Ba r ) < 2 are formed. 

We conclude that NSM have to respect stringent re- 
quirements if they are to operate as dominating r-process 
contributors, which most likely are not fulfilled in reality: 
High NSM rates of the order > 10 _3 yr _1 are required 
to account for the early appearance of r-process nuclei 
in the ISM. Consequently, the amount of r-process matter 
ejected in a NSM event has to be of the order < 1CP 4 Mq, 
lest the total inventory of r-process nuclei in the Galaxy 
is exceeded. In addition, the coalescence timescale of a 
large fraction of neutron star binaries has to be of the or- 
der 1-10 Myr or at the most 100 Myr. If these constraints 
are indeed fulfilled in reality, then the properties of NSM 
(i.e. rates and ejected r-process matter) are comparable to 
the ones of SNe II that (may) give rise to r-process nucle- 
osynthesis from the point of view of chemical evolution. 
Consequently, the enrichment of the ISM with r-process 
matter from NSM would be similar to the one where core- 
collapse SNe are the source of r-process elements (compare 
Figs. n~Elwhh the upper left panels of Figs. 0] and [5J. 

However, it seems unlikely that NSM fulfil these con- 
ditions, which suggests that NSM are ruled out as the 
major r-process source. On the other hand, NSM occur- 
ring at low rates and with low ejecta masses of the order 
10 -4 — 10~ 5 Mq may still contribute to the enrichment of 
the ISM with r-process nuclei if their r-process signature 
is similar to the one generated by SNe II. In this case, 
the impact of NSM on r-process nuclei enrichment would 
be negligible compared to the fast injection of r-process 
matter by SNe II. 

4. Conclusions 

In this work, we study the enrichment of the interstel- 
lar medium with r-process elements in the framework of 
inhomogeneous chemical evolution. We present a detailed 
comparison of the impact of lower-mass SNe II (8 — 10 Mq, 
model SN810), higher-mass SNe II (> 20 Mq, models 
SN2025 and SN2050) and NSM as major r-process sites 
on the enrichment history of the early Galaxy. 

Among the SNe II scenarios, model SN2050 gives the 
best fit to observations, since it reproduces the trend of 
r-process abundances at ultra low metallicities and at the 
same time fits observations with the highest [r/Fe] abun- 
dance ratios at [Fe/H] « —3. Neither model SN2025 nor 
SN810 are able to reproduce observations of the most r- 
process enriched stars. In addition, model SN810 gives the 
worst fit to the trend of r-process abundances at ultra low 
metallicities. However, we conclude that, due to the large 



uncertainties inherent in the progenitor mass dependence 
of iron yields of SNe II, it is not possible to clearly rule out 
either lower-mass SNe II or higher-mass SNe II as dom- 
inant r-process sites from the point of view of inhomo- 
geneous chemical evolution. Additional uncertainties are 
introduced by the fact that reliable r-process yields from 
SNe II are unavailable as yet. Here, they were deduced in 
such a way that the average [r/Fe] abundances in metal- 
poor halo stars are reproduced. Clearly, iron and r-process 
yields and their dependence on progenitor mass from self- 
consistent core-collapse SN calculations are required be- 
fore any decisive conclusion can be reached. 

On the other hand, NSM seem to be ruled out as major 
r-process sources for the following reasons: 

1. Estimates of the galactic NSM rate are two to three 
orders of magnitude lower than estimates of galactic 
SNe II rate. Thus, the injection of r-process nuclei into 
the interstellar medium by NSM would occur very late 
during Galaxy formation ([Fe/H] w —2.5), whereas r- 
process elements are already observed at [Fe/H] = —3 
and probably even at [Fe/H] = —4. 

2. The late injection of r-process elements furthermore 
leads to prominent tails in the distribution of r-process 
abundances down to very low [r/Fe] ratios at [Fe/H] 
> — 2, which are seemingly not consistent with obser- 
vations. It is possible that r-process abundances [r/Fe] 
< —1 are not detectable to date in metal-poor halo 
stars. However, the large fraction of stars predicted by 
our model with -1 < [r/Fe] < at [Fe/H] > -2 should 
be present in the observational data. 

3. Since NSM occur at a lower rate than SNe II, their 
r-process yield has to be about two orders of magni- 
tudes higher than the r-process yield of typical SNe II. 
Due to this large r-process yield, considerable chemical 
inhomogencities in the ISM are expected to be present 
even at solar metallicity. The scatter in [r/Fe] is pre- 
dicted to be of the order 2.0—2.5 dex, whereas a scatter 
of only 0.2 — 0.3 dex is observed. 

NSM as major r-process sources are only consistent 
with observations under the following conditions: First, 
they occur at rates about one to two orders of magni- 
tude higher than given by present estimates. Second, a 
dominant part of the NSM population has coalescence 
timescales shorter than approximately 10 Myr. Third, the 
total amount of ejected r-process matter is of the same 
order as the present theoretical estimate from relativistic 
merger calculations. Under these conditions, the enrich- 
ment of the ISM with r-process elements from NSM is 
qualitatively very similar to the enrichment of r-process 
elements dominated by SNe II. However, it seems highly 
unlikely that NSM fulfil these conditions, which suggests 
that NSM are ruled out as the major r-process source. 
Nevertheless, it is still feasible from the point of view of 
chemical evolution that NSM contribute a minor fraction 
of r-process elements to the enrichment of the ISM. In this 
case they have to occur at a low rate and the ejected r- 
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process matter must not exceed the lower limit indicated 
by recent relativistic merger calculations. 

Thus, we conclude, that the exact astrophysical nature 
of r-process sites still remains a mystery, since it is not 
possible to clearly distinguish between neutron capture 
element abundance patterns resulting from lower-mass 
SNe II (8 - 10 M ) and the ones from higher- mass SNe II 
(> 20 Mq) in the framework of inhomogeneous chemi- 
cal evolution. However, the present investigation suggests 
that core-collapse SNe are much more likely to be the 
dominant r-process sites than coalescing neutron star bi- 
naries, which at least reduces the list of possible major 
contributors of r-process nuclei to the enrichment of the 
interstellar medium. Yet, it remains to be seen how SNe II 
can actually produce the required r-process yields. 
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